%To compare the transport  estimates from Voltage measurements
% across florida state with the ECCO electric field simulations
% Used data from The florida cable transport project at
% aoml.noaa.gov/phod/floridacurrent/index.php
%latest date 6 Nov 2005

load c:\manoj\ocean\velocities\UVECCO.mat IV TimeID; %loading the IV -> as the floe across florda is pred NS
load c:\manoj\ocean\velocities\florida\transport fld datefl fluv dateflorida; %fld -> transport estimated florida(:,5) is the measured transport
load c:\manoj\ocean\x3dg\cable_results\5Nov05_old_fmt.mat Etr Epr; % The simulated Electric field from ECCO 1992-2002 fields

%Notes
%IU is irrelevant at least 1 times smaller than IV in this region
%FLORIDA STRAIT
clat = [26.5,26.5]; % start and end lat for florida measurements (only guess !)
clon = [-75.7,-80.0];

%DRAKE PASS

% worldmap(flipud(squeeze(IU(1,:,:))),[1,90,-180]);
% worldmap(flipud(squeeze(Etr(1,:,:))),[1,90,0]); % The last 0 was checked several times !

dateE = TimeID/(3600*24) + datenum(1992,1,0);

for i = 1:132,
[V,r1,lat,lon] = mapprofile(flipud(squeeze(Epr(i,:,:)))*1e6,[1,90,0],clat,clon);
Vol(i) = trapz(deg2km(r1),V);
end;

for i = 1:132,
[V,r1,lat,lon] = mapprofile(flipud(squeeze(IU(i,:,:))),[1,90,-180],clat,clon);
Sv(i) = trapz(deg2km(r1)*1e3,V);
end;

plot(dateflorida,fluv(:,5),'b.');
hold on;
plot(datefl,fld(:,4),'b');
plot(dateE,Sv/1e6,'r-');

% Note the ve 